Stable rotating dipole solitons in nonlocal optical media 
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We reveal that nonlocality can provide a simple physical mechanism for stabi- 
lization of multi-hump optical solitons, and present the first example of stable 
rotating dipole solitons and soliton spiraling, known to be unstable in all types 
of realistic nonlinear media with local response. 
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Recently increased interest in the study of self-trapped 
optical beams in nonlocal nonlinear media is explained 
by experimental observations of nonlocal spatial solitons 
in liquid crystals^ and lead glasses^ as well as a num- 
ber of many interesting theoretical predictions^i^, in- 
cluding stabilization of vortex solitons against symmetry 
breaking instability in the media with nonlocal optical re- 
sponse 7 '^. Many of the predicted and demonstrated prop- 
erties of nonlocal nonlinear modes suggest that in such 
optical media we should expect stabilization of many dif- 
ferent types of spatial nonlinear structures such as neck- 
laces^, soliton clusters^, and recently introduced broader 
classes of modulated optical vortices or azimuthonsii. 

One of the simplest multi-hump solitons predicted 
in nonlinear optics is a dipole-like structure composed 
of two interacting fundamental beams with opposite 
phases that undergo angular rotation during the prop- 
agationiSiiii. Such a localized structure can be viewed as 
a special type of strongly modulated single-charge vor- 
tex soliton or azimuthon and it can also be linked to the 
problem of soliton spiraling, as discussed e.g. in §6.1 of 
Ref. [13]. However, in spite of its broad applicability and 
deep physical context, the rotating dipoles are known to 
be always unstable in all types of realistic optical media 
with local respons o 12 i 13 . 

In this Letter, we study the effect of nonlocal response 
of an optical medium on the stability of rotating dipole 
solitons. We demonstrate that such solitons can be sta- 
bilized when the nonlocality parameter exceeds some 
threshold value that decreases with the angular momen- 
tum and angular velocity, i.e. the dipoles rotating faster 
are more robust. We describe a continuous family of two- 
lobe localized structures (or azimuthonsii) that extend 
from non-rotating dipole solitons to radially symmetric 
vortex solitons by continuous azimuthal transformations. 

We consider the propagation of paraxial optical beams 



in a nonlinear medium described by the nonlinear 
Schrodinger equation for the scalar field envelope £— 
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n(I,r)E = 0, 



(1) 



where z and r stand for the propagation and transverse 
coordinates, respectively. We assume that the nonlinear 



refractive index n depends on the intensity / 
the following nonlocal relation 



n(I, r) 



R(\r-p\)I(p)dp, 



\Er via 
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where the response function R(r) is determined by the 
specific physical process responsible for the medium non- 
linearity. Here we consider a Gaussian nonlocal response 



R(r) = a^n- 1 exp (-r 2 /a 2 ) , 



(3) 



where a measures the degree of nonlocality. In the fol- 
lowing we use the scaling f= fa, z — z'a 2 , E = E'/a, 
and omit primes. 

Multi-hump localized solutions of this model can be 
found numerically by two-dimensional relaxation meth- 
odsi*. However, this approach represents a major dif- 
ficulty for the complex envelopes with singular phase 
profiles due to a poor convergence. Therefore, here we 
employ an approximate but much simpler variational 
method 15 allowing a reasonably good approximation 11 . 
To this end, we introduce the following ansatz 
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^ttU (r) [cos (my;) + ip sm(rrvp)] exp (ikz) , (4) 



where ip is the azimuthal angle and the parameter p (0 < 
p < 1) determines the modulation depth (contrast) of 
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Fig. 1. (Color online) Beam radial intensity U 2 (r) (solid) 
and effective nonlinear potential N(U, r) (dashed) for the 
dipole soliton (p = 0, red) and vortex soliton (p = 1, 
blue) for (a) relatively low nonlocality with k = 1 and 
(b) high nonlocality with k = 70. (c) Power vs. prop- 
agation constant k for the dipole solitons; several lines 
for different p coincide, (d)-(f) Envelopes of the intensity 
(top) and phase (bottom) for three characteristic cases 
with k = 70 and (d) p = 0, (e) p = 2" 1 / 2 , and (f) p = 1. 

the beam intensity, n = 1 — p 2 . Substituting this ansatz 
into the averaged Lagrangian 15 associated with Eq. Q), 
we obtain the equation for the radial amplitude U(r) of 
the localized mode, 

-W + U rr + r~ l U r - m 2 r- 2 U + UN{U 2 ,r) = 0, (5) 

where the nonlinear term is given by the expression 

/>oo 

N(U 2 ,r) = ir{l+p 2 )c- r2 / pdpe- p2 U 2 (p) x 

Jo 

Io(2rp) + 1 (j^j hm(2rp)\ ,(6) 

and Ij is the modified Bessel function of the first kind. 
In the limit p — > 1, the function (J2J becomes radially 
symmetric, and the potential © corresponds to a vortex 
soliton. On the other hand, for p — > 0, this solution 
describes scalar multipole solitons (or optical necklaces) 
known to be radially unstable in local media. 

Applying to Eq. J^J) the numerical shooting technique 
and an iteration procedure, we find the stationary local- 
ized solutions with different values of p and k. Our results 
are summarized in Figs.^a-f). Figures^a,b) show the 
radial envelopes of the localized modes (solid lines) and 
the corresponding profiles of the effective refractive index 
(dashed lines) for the dipole (p = 0) and vortex (p = 1) 
solitons in the case of low [(a)] and high [(b)] degree of 
nonlocality. Figure^c) shows the dependence of the soli- 
ton power P = tt(1 + p 2 ) f °° U 2 rdr on the propagation 
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Fig. 2. (Color online) Two different instability scenarios 
for rotating dipole solitons with p = 0.5 for (a) k = 1 
and (b) k = 20. (c) Stable rotation for p — 0.5 and 
k = 100. In all three cases initial noise is 10%. (d) 
Stability domain of dipole solitons in the (p, k) space, (e) 
Numerically found relation between the angular velocity 
Q of a dipole soliton and its contrast parameter p. 



constant fc, for several values of the modulation param- 
eter p. Figures ^d-f) present the intensity (top) and 
phase (bottom) distributions for the dipole, azimuthon, 
and vortex solitons, respectively. 

Having a nontrivial phase structure, the azimuthons 
carry a nonzero angular momentum. Using the ansatz 
Eq. 0, we find that the beam orbital angular momen- 
tum, M = Im J E*d v E dr, normalized to its power P, 
defines the beam effective (fractional) spin, S = M/P = 
2p/(l+p 2 ). 

These results show the effect of the modulation depth 
(determined by p) on the soliton intensity, and they 
demonstrate that a transition from a dipole mode to a 
radially symmetric vortex soliton can be achieved by con- 
tinuous azimuthal transformations. Therefore, the two- 
peak azimuthon solutions described above provide a link 
between two spiraling solitons with opposite phase and 
radially symmetric optical vortices, or spatially localized 
singular beams. 

In order to study stability of the stationary solutions 
we solve Eq. Q numerically by the split-step beam prop- 
agation method with a fast Fourier transform. As initial 
conditions, we take the stationary solutions perturbed 
initially by 10 — 20% of noise and propagate them over 
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Fig. 3. (Color online) (a,b) Evolution of a dipole soliton 
(k = 60, p = 0) excited by a strong initial deformation. 
Persistent oscillations of both the mode amplitude / max 
(dashed) and diameter D (a distance between two peaks, 
solid) are observed; these oscillations decay slow during 
the propagation, (c) Example of the soliton revivals ob- 
served in the instability-induced evolution of double-ring 
single-charge vortex solitons with k = 100. 



modulational instability is fully developed, as shown in 
Fig. 3(c). Such an instability can be easily excited if we 
perturb the stationary localized mode by a small amount 
S as follows, U(r)[cos(mip) + ism(rrvp + 6)]. The struc- 
tures which appear at intermediate stages of this insta- 
bility dynamics are similar to higher-order linear modes 
such as Laguerre- and Hermite-Gaussian modes. Thus, 
in this case the nonlocal medium acts as a self-induced 
mode converter, where the conversion dynamics depends 
on the initial mode and the properties of the medium, 
i.e. the value of the nonlocality parameter. After several 
revivals, the structure losses its symmetry and undergoes 
irregular dynamics. 

In conclusion, we have demonstrated that nonlinear 
optical media with a nonlocal response can stabilize ro- 
tating dipole solitons and their generalizations in the 
form of modulated azimuthons with a fractional spin. 
The stabilization is achieved when the nonlocality param- 
eter exceeds a certain threshold value, and faster rotating 
azimuthons are found to be more robust. 

S. Lopez- Aguayo is with Photonics and Mathematical 
Optics Group, Tecnologico de Monterrey, Mexico. 



the distances of z ~ 10 2 — 10 3 or larger. 

Figure El presents typical propagation scenarios. When 
the nonlocality is weak the dipole breaks into two mutu- 
ally repelling filaments [(a)], and this behavior is rem- 
iniscent of the azimuthal instability of vortices in local 
mediai 3 -. For the subcritical nonlocality [(b)] the az- 
imuthon initially breaks up into two filaments which are 
subsequently forced by the nonlocality-induced potential 
to collide and merge into a single fundamental soliton. 
Stable propagation of the rotating azimuthon is observed 
when nonlocality exceeds a critical value [(c)] . This corre- 
sponds to a balance of the effective nonlocality-mediated 
attractive force and centrifugal repulsive force. The do- 
main of stability on the plane (p, k) is shown in Fig.[5fd), 
and the stability threshold decreases in the limit of the 
radially symmetric vortices (p — > 1). For the stable ro- 
tating dipoles, we obtain the angular velocity by nu- 
merical averaging and Fig. Hfe) shows the dependence of 
the azimuthon angular velocity on the modulation pa- 
rameter p. Velocity increases monotonically with 5* and 
only weakly depends on the power (two curves are for 
k = 60 and k = 80). For the vortex soliton with p = 1 
the rotation velocity is a free parameter. 

We have tested the stability of the rotating dipole 
solitons with respect to strong perturbations which lead 
to large-amplitude oscillations, as shown in Fig. 3(a,b). 
Similar studies have been carried out for the spatial soli- 
tons described by Eq. JIJ with higher azimuthal indices 
m > 1 (not shown). Somewhat surprising results we 
obtained by testing stability of double-ring radially sym- 
metric (p = 1) solitons with m = 1. While no absolutely 
stable solitons have been found in our numerical simula- 
tions, we have observed the effect of soliton revivals after 
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